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The reduced cross sections for ep deep inelastic scattering have been measured 
with the ZEUS detector at HERA at three different centre-of-mass energies, 318, 
251 and 225 GeV. From the cross sections, measured double differentially in 
Bjorken x and the virtuality, Q^, the proton structure functions and F2 have 
been extracted in the region 5 x 10"^ < x < 0.007 and 20 < < 130 GeV^. 
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1 Introduction 



The inclusive e^p deep inelastic scattering (DIS) cross section can, at low virtuality of 
the exchanged boson, Q^, be expressed in terms of the two structure functions, F2 and 

Fl, as 

= (^ix,Q ,y), (1) 

where a is the fine structure constant, x is the Bjorken scaling variable, y is the inelasticity 
and = 1 + (1 — y)^. The quantity a is referred to as the reduced cross section. The 
kinematic variables are related via = xys, where ^/s is the ep centre-of-mass energy. 
The magnitude of is proportional to the absorption cross section of longitudinally 
polarised virtual photons by protons, Fl oc ai, while F2 includes also the absorption 
cross section for transversely polarised virtual photons, F2 oc {ar + o'l)- At low values of 
X and small o"/,, the ratio R = Fl/{F2 — F^) ^ ctl/o't gives the relative strength of the 
two components. 

HERA measurements of the reduced ep DIS cross section and F2 [1-3] provide the 
strongest constraints on the proton parton distribution functions (PDFs) at low x. Within 
the DGLAP formalism [4], F2 at low x is dominated by the qq sea distributions while the 
scaling violations of F2 reflect the gluon distribution, g{x, Q'^), via a convolution with the 
splitting function Pqg{x), dF2/d\ia.Q^ ~ as{Q^)Pqg{x) <^xg{x,Q^), where as is the strong 
coupling constant. 

The published values of F2 at low x at HERA required assumptions to be made about Fl 
or were restricted to the kinematic region where the contribution from F^ was sufficiently 
suppressed to be neglected. Moreover, gluon distributions extracted from scaling viola- 
tions are dependent on the formalism [5] and the order of perturbative expansion [6] used 
to calculate the splitting functions. Measurements of the reduced cross section at fixed 
{x, Q"^) and different y allow F2 and Fl to be extracted simultaneously, thereby eliminating 
the assumptions about Fl when extracting F2. Furthermore, a direct measurement of Fl, 
which is strongly correlated to the gluon density [7], provides an important consistency 
check of the formalism. 

A model-independent determination of Fl requires the reduced cross section to be mea- 
sured at fixed values of x and for multiple centre-of-mass energies (varying y values). 
This method has been previously used to extract Fl in fixed-target experiments [8-11] and 
recently by the HI collaboration [12]. The HI collaboration has also applied extrapolation 
methods to determine Fl [2,13]. 

In this paper, the first ZEUS measurement of Fl is presented as well as the most precise 
ZEUS measurement of F2 in the kinematic region studied. Comparisons of theoretical 



dxdQ'^ 



xQ' 



F2{x,Q')-^Fl{x,Q') 

-I 4~ 
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predictions with the data are also presented. 



2 Experimental method 

The values of F2 and were extracted at fixed x and by fitting a straight line 
to the values of a against jY^ in the so-called Rosenbluth plot [14]. The method 
is based on Eq. [1], which implies that F^ix.Q'^^ = a{x,Q'^,y = 0) and Fl{x,Q'^) = 
—da{x,Q'^,y)/d{y'^/Y^), which in turn implies the need for data at fixed (x, Q^) and 
different y. At HERA, this can be achieved by varying y/s. 

The precision of this procedure depends on the range available in y'^/Y^. This was 
maximised by collecting data at the nominal HERA energy, ^/s = 318 GeV, and at 
a/s = 225 GeV, the lowest attainable energy with adequate instantaneous luminosity. 
An intermediate data set was collected at y/s = 251 GeV. 

The variation of ^/s was achieved by varying the proton beam energy, Ep, while keeping 
the electron beam energy constant, E^. = 27.5 GeV. Data were collected in 2006 and 
2007 with Ep = 920, 575 and 460 GeV, referred to respectively as the high- (HER), 
medium- (MER) and low-energy-running (LER) samples. The corresponding integrated 
luminosities of the HER, MER and LER samples are 44.5, 7.1 and 13.9 pb~^, respectively. 

3 Experimental apparatus 

A detailed description of the ZEUS detector can be found elsewhere [15]. A brief outline 
of the components most relevant for this analysis is given below. 

In the kinematic range of the analysis, charged particles were tracked in the central 
tracking detector (CTD) [16] and the microvertex detector (MVD) [17]. These components 
operated in a magnetic field of 1.43 T provided by a thin superconducting solenoid. The 
CTD drift chamber, consisting of 72 sense wire layers organised into 9 super layers, covered 
the polar- angl^ region 15° < 6* < 164°. The MVD silicon tracker consisted of a barrel 
(BMVD) and a forward (FMVD) section. The BMVD provided polar-angle coverage for 
tracks with three measurements from 30° to 150°. The FMVD extended the polar-angle 
coverage in the forward region to 7°. 

The high-resolution uranium-scintillator calorimeter (CAL) [18] consisted of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part was 

^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing towards the 
centre of HERA. The coordinate origin is at the nominal interaction point. 
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subdivided transversely into towers and longitudinally into one electromagnetic section 
(EMC) and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections (HAC). 
The smallest subdivision of the calorimeter was called a cell. The CAL energy resolutions, 
as measured under test-beam conditions, were a{E)/E = 0.18/ \/E for electrons and 
(r{E)/E = 0.35//E for hadrons, with E in GeV. 

The rear hadron-electron separator (RHES) [19] consisted of a layer of approximately 
10 000 (3x3 cm^) silicon-pad detectors inserted in the RCAL at a depth of approximately 
3 radiation lengths. The polar-angle coverage is approximately 131° < 9 < 173°. The 
small-angle rear tracking detector (SRTD) [20] was attached to the front face of the 
RCAL and consisted of two square planes of scintillator strips. The detector covers the 
total area of 68 cmx68 cm, with a 20 cmx20 cm cutout in the centre for the beam-pipe. 
The polar-angle coverage is 162° < 9 < 176°, with full acceptance for 167° < 9 < 174.5°. 

The small tungsten-scintillator calorimeter located approximately 6 m from the interac- 
tion point in the rear direction was referred to as the "6m-tagger" [21]. For scattered 
electrons in the energy range from 4.1 to 7.4 GeV, the acceptance was very close to one 
with very high purity. 

The luminosity was measured using the Bethe-Heitler reaction ep — > e^p by a luminosity 
detector which consisted of an independent lead-scintillator calorimeter [22] and a mag- 
netic spectrometer [23] system. The fractional systematic uncertainty on the measured 
luminosity was 2.6%. 



4 Event reconstruction and selection 

The kinematic region studied spanned 0.09 < y < 0.78 and 20 < < 130 GeV^, corre- 
sponding to 5 X 10~^ < X < 0.007. The event kinematics were evaluated based on the 
reconstruction of the scattered electron [24] using 

2/e= 1-^(1- COS^e), (2) 

Ql = 2ElE,{l + cos9e), (3) 

where 9^ and E'^ are the polar angle and energy of the scattered electron, respectively. 

Electrons were identified using a neural network based on the moments of the three- 
dimensional shower profile of clusters found in the CAL [25]. The quantity E'^ was re- 
constructed using the CAL, and 9e was determined using the reconstructed interaction 
vertex and scattered electron position in the SRTD or, if outside the SRTD acceptance, 
in the RHES. In less than 2% of events, 9e could not reliably be determined using the 
SRTD+HES system and such events were rejected. 
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The quantity 6 = ^^{E — pz)i was used both in the trigger and in the offline analysis. 
The sum runs over all CAL energy deposits. Conservation of energy, E, and longitudinal 
momentum, pz, implies that 6 = 2Ee = 55 GeV. Undetected particles that escape through 
the forward beam-pipe hole contribute negligibly to 6. Undetected particles that escape 
through the rear beam-pipe hole, such as the final-state electron in a photoproduction 
event, cause a substantial reduction in 6. Events not originating from ep collisions often 
exhibit a very large S. 

A three-level trigger system was used to select events online [15,26-28]. A dedicated 
trigger was developed providing high efficiency for high-y events [29]. The trigger required 
an event to have 6 > 30 GeV and either an electron candidate with E'^ > 4 GeV in the 
RCAL and outside a 30 cmx30 cm square centred around the beam-pipe, or ^^<i65° ^ 
20 GeV, where ^^<i65° (ienotes 6 calculated only from the CAL energy deposits at polar 
angles less than 165°. 

Events were selected offiine if: 

• 42 < 5 < 65 GeV; 

• the reconstructed interaction vertex fulfilled |^vtx| < 30 cm; 

• the energy of the most probable electron candidate satisfied E'^ > Q GeV; 

• the event topology was not compatible with an elastic QED Compton (QEDC) event; 

• the event timing was consistent with the HERA bunch structure; 

• He < 0.95 and yjB > 0.05, where i/jb is the Jacquet-Blondel estimator [30] of y; 

• PT,h/pT,e > 0.3, where pT,h and pT,e refer to the transverse momentum of the hadronic 
system and electron candidate, respectively. 

The projected path of the electron candidate was required to: 

• exit the CTD at a radius > 20 cm and hence traverse the MVD fiducial volume and 
at least four CTD sense-wire layers, ensuring the possibility of identifying the track; 

• enter the RCAL at a radius < 135 cm, missing the inactive region between the RCAL 
and BCAL sections. 

Hit information from the MVD and CTD was used to identify the tracks of the elec- 
tron candidates. The procedure was based on the ratios of the number of observed to 
the maximum number of possible hits in the MVD and CTD, denoted f^^^ and /hi™, 
respectively: 



?MVD 
hit 



> 0.45; 



?CTD 
hit 



> 0.6. 
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This method was used because of the wider polar angular acceptance compared to the 
regular tracking capability of the MVD+CTD tracking system. Specifically, for an event 
with a nominally placed vertex, the electron candidate can be validated up to an angle of 
9e = 169 o, compared to 6e = 159 o with full tracking. After all cuts the HER, MER and 
LER samples contained 819168, 115719 and 205967 events, respectively. 



5 Cross section determination 

The reduced cross sections in a given (x, Q^) bin were calculated according to 

AT _ AT^^ 

) = ttdIs — ^^SM(a;,Q ), 

where 5"sm(3^, is the Standard Model electroweak Born- level reduced cross section and 
A^data, -^MC ^^"^ -^MC denote, respectively, the number of observed events in the data and 
the expected number of background and DIS events from the Monte Carlo (MC). The 
CTEQ5D [31] parameterisation of the proton PDF was used when calculating a"sM(2:, Q^) 
as well as in the MC models when evaluating and N'^q- Specifically, the DIS signal 
processes were simulated using the Djangoh 1.6 [32] MC model. After the full event 
selection, the background consisted almost entirely of photoproduction events. These 
were simulated using the Pythia 6.221 [33] MC model. The additional background 
components that were considered were elastic QEDC and mis-reconstructed low-Q^ DIS, 
simulated using the Grape-Compton [34] and Djangoh 1.6 MC models, respectively. 
The MC events were processed through a full simulation of the ZEUS detector and trigger 
system based on Geant 3.21 [35]. 

The Djangoh and Pythia samples included a diffractive component and first-order 
electroweak corrections. The diffractive and non-diffractive components of the Djangoh 
sample were scaled to improve the description of the HER, MER and LER ?7niax distribu- 
tions, where ?7niax is equal to the pseudorapidity of the most forward CAL energy deposit. 
The electroweak corrections were simulated using the Heracles 4.6 [32,36] MC model. 
Their uncertainty was evaluated by comparing the predictions from Heracles to the 
higher-order predictions from HECTOR 1.0 [37]. The predictions were found to agree 
to within 0.5%. The hadronic final state of the Djangoh MC was simulated using the 
colour-dipole model of Ariadne 4.12 [38] which uses the Lund string model of Jetset 
7.4 [39] for the hadronisation. 

In order to improve the Monte Carlo description of the photoproduction component, the 
contribution from the direct subprocesses was enlarged from 3% (default) to 9% in the 
inclusive Pythia sample while contributions by diffractive subprocesses were reduced 
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accordingly. This procedure ensured that previous ZEUS resuhs were reproduced [40,41] 
and the predicted inclusive Pythia cross section remained unchanged. The predicted 
photoproduction cross sections for HER, MER and LER were then validated against 
photoproduction data samples selected using the 6m-tagger. The predicted cross sections 
were consistent with these data within the ±10% total uncertainty on the data. 

Figure [1] shows the distributions of the variables E'^ and 6e within the HER, MER and 
LER data sets compared to the combined detector-level predictions from the MC models. 
The agreement is good in all cases. According to the MC models, the final data sample 
contained 97% DIS signal and 3% background events. The vast majority of the background 
events were found at low and high y; in the most affected kinematic bin, the background 
fraction was 16%. 

The reduced cross sections, a, were measured from the HER, MER and LER samples in 
the kinematic region 0.09 < y < 0.78 and 20 < < 130 GeV^. The a are given double 
differentially in x and in Tables [TH3l The a are also shown at the 6 selected values 
as functions of x in Fig. [21 The cross sections have been compared to DGLAP-predictions 
based on the NLO (C(a^)) ZEUS- JETS PDF set [42], as well as the prediction for Fl = 0. 
The QCD prediction with a non-zero F^, describes the data well and is favoured over 
Fl = 0. 

6 Systematic uncertainties 

The systematic uncertainty on the reduced cross sections due to the following sources were 
evaluated [29] (the numbers in the parentheses are the maximum uncertainty observed in 
any one of the reduced cross section bins): 

• {^-/p}^ ±10% uncertainty on the level of photoproduction background (—2%); 

• {^Ee}: the electron energy-scale uncertainty of ±0.5% for E'^ > 20GeV, increasing to 
±1.9% at E'^ = 6GeV(4.4%); 

• {Se,,}, the ±2% hadronic energy-scale uncertainty (—4.1%); 

• {^cId}, the uncertainty on the electron-finding efficiency, evaluated by loosening (tight- 
ening) the criterion applied to the output of the neural network used to select electron 
candidates, both for data and MC (±1.8%); 

• {^dx, ^dy}, the SRTD and HES position uncertainty of ±2 mm in both the horizontal 
and vertical directions (±3%); 

• {<^MVD5 ^ctd}, the uncertainty on the hit-finding efficiency, evaluated by loosening 
(tightening) the hit fraction criteria, both for data and MC (±3.7%); 
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• {(^diff}, the ±10% uncertainty on the scale factors apphed to the diffractive Djangoh 
component (—0.7%). 

The one-standard-deviation systematic uncertainties due to each source are hsted in 
Tables [TH3] for the reduced cross sections at the three different centre-of-mass energies. 
All of the uncertainties are symmetric. They are quoted with a sign indicating how the 
reduced cross sections would vary given an upwards variation in the electron or hadronic 
energy scales, the SRTD and HES positions, the photoproduction cross section or the 
diffractive scale factors, or looser selection criteria on the neural network output or MVD 
or CTD hit fractions. 

The total uncertainty on the normalisation included 

• the luminosity uncertainty, which was ±2.6% for all three data sets, of which ±1% 
was uncorrelated between the data sets; 

• the uncertainty on simulating the interaction-vertex distribution, evaluated by com- 
paring the ratio of the number of events with {Z^txl < 30 cm and {Z^txl > 30 cm in 
data and MC (±0.3%); 

• the trigger-efficiency uncertainty (±0.5%). 

The luminosity, vertex-distribution and trigger-efficiency uncertainties are perfectly corre- 
lated between bins and hence, when added in quadrature, constitute a total normalisation 
uncertainty of ±2.7%, of which ±2.5% was correlated between the running periods and 
±1.1% uncorrelated. The uncertainty due the electroweak corrections was found to be 
negligible. 

The total systematic uncertainty in each bin, formed by adding the individual uncer- 
tainties in quadrature, is also given in Tables [IH3l This sum also includes the statistical 
uncertainty due to the combined MC sample {Sunc} and is the only systematic source that 
is considered to be uncorrelated between bins. This total systematic uncertainty does not 
include the total normalisation uncertainty. Propagation of the systematic uncertainties 
to Fl, F2 and R is described in the next section. 

7 Extraction of Fl, F2 and R 

In order to extract Fl , F2 and R a different binning scheme than that given in Tables [THU 
was applied to the reduced cross sections. Bins in y were chosen such that, for each of the 
6 bins, there were 3 values of x at which the reduced cross sections were measured from 
all three data sets. This removed the need to interpolate the data between different points 
in the [x, Q"^) plane. The structure functions were extracted by performing a simultaneous 
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fit to these 54 measured cross section values using Eq.[Tl Prior to fitting, the three data sets 
were normahsed to their luminosity-weighted average in the restricted kinematic region, 
y < 0.3, where the contribution to the reduced cross sections from is small. This 
procedure resulted in scaling the data by factors of 1.0027 ± 0.0027, 0.9869 ± 0.0051 and 
0.9997 ±0.0039, for the HER, MER and LER data sets, respectively. The spread of these 
factors is consistent with the uncorrelated part of the total normalisation uncertainty of 
1.1%. 

To extract Fl and F2, 48 parameters were fit simultaneously: 18 F2 and 18 Fl values for 
the 18 (x, Q^) points; 3 relative normalisation factors for the HER, MER and LER data 
sets and 9 global shifts of systematic uncertainties (5-yp, 5e^-, ^e^i ^eiD, ^dx, ^dy, ^mvd, ^ctdi 
(5diff)- The three normalisation factors allowed for variations of the relative normalisation 
factors within their remaining uncertainties (see above). The nine global shifts allowed for 
changes in the central values of a in a correlated manner across the (x, Q'^) plane according 
to the uncertainties listed. The probability distributions for the shifts of the systematic 
sources and the relative normalisations were taken to be Gaussian, with standard devi- 
ations equal to the corresponding systematic uncertainty. The probability distributions 
for the cross sections at each (x, Q"^) point were also taken to be Gaussian with standard 
deviations given by 5stat and 5unc added in quadrature. The fit was performed within the 
BAT (Bayesian Analysis Toolkit) package [43] which, using a Markov chain MC, scans 
the full posterior probability density function in the 48-dimensional parameter space. 

Initially, the F^ and F2 parameters were left unconstrained and fiat prior probabilities 
were assumed. The results are given in Table HJ and are labelled with the superscript (1). 
The values quoted in the table were evaluated at the point where the probability den- 
sity function attains its global maximum. The uncertainty ranges correspond to minimal 
68% probability intervals. These ranges represent the full experimental uncertainty, which 
comprises statistical as well as systematic uncertainties. The fitted shifts, representing the 
correlated variation of the data points according to relative normalisation and correlated 
systematic uncertainties, are typically within 0.1 and at most 0.5 standard deviations of 
the normalisation or systematic uncertainties. The F2 values typically have uncertainties 
of 0.03, while the Fl values have uncertainties ranging from 0.1 to 0.2. These F2 mea- 
surements are the most precise available from the ZEUS collaboration in the kinematic 
region studied here. The results are shown in Fig. [3] together with predictions from the 
ZEUS- JETS PDF fit. Good agreement is observed. 

Applying constrained priors -F2 > and < F/, < F2 in the fitting gave marginally 
different results as seen in Table H] (results are denoted with the superscript (2)). For 
example, the most probable value for Fl at = 45GeV^ and x = 0.00153 is now 0, in 
which 68% probability upper limit is given. 

Further fits to the data were performed to extract Fl{Q'^)i R{Q^), and a single overall 
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value of R for the full data set. In each case, the same fitting procedure as described 
above was used, but with a reduced number of parameters. To extract Fl(Q^), first 
r{Q'^) was fitted, where r = Fi^jF^- In fitting r[Q'^\ a single value of r was taken for 
all X points in the same bin. Only a weak dependence of r on x in a restricted x 
range is expected in the NLO DGLAP formalism as well as in phenomenological models. 
This prediction is supported by the data within the experimental uncertainties. Flat prior 
distributions for r((5^) were assumed and both unconstrained and constrained fits were 
made, with r(Q^) > enforced in the latter. The value of Fl{Q'^) was then evaluated 
as Fi^Xi^Q"^) = r{Q'^)F2{xi,Q'^), where for each point, Xj was chosen such that Q'^/xi 
was constant, which for a/s = 225 GeV, corresponds to y = 0.71. The results are given 
in Table [5] and the unconstrained values are shown in Fig. H^. These data are in good 
agreement with the results obtained by the HI collaboration [12]. 

Values of R{Q^) and an overall value of R were extracted with fiat prior distributions. 
Both unconstrained and constrained fits were made. In the latter, it was required that 
< Fl{Q^) < F2{Q^) and < < Fg. The results from both fits are given in Table [6] 
and the unconstrained R{Q'^) values are shown in Fig. |3)d. The uncertainty in the overall 
R is not reduced as much as might be expected compared to the uncertainties on R{Q^) 
due to the correlation between the values at different . The value of R from both the 
unconstrained and constrained fits was R = O.ISIqqj. 

Figures and also show a comparison of the data with predictions based on the 
ZEUS- JETS and CTEQ6.6 [44] NLO and MSTW08 [45] NLO and NNLo|fits. All these 
predictions are based on the DGLAP formalism^]. Also shown are predictions from the 
NLL BFKL resummation fit from Thorne and White (TW) [48] , and the prediction from 
the impact-parameter-dependent dipole saturation model (b-Sat) of Kowalski and Watt 
based on DGLAP evolution of the gluon density [49]. All of the models are consistent 
with the data. 



8 Summary 

The first measurement of Fl{x, Q^) by the ZEUS collaboration is presented, as is the first 
measurement of F2{x, Q"^) at low x that does not include assumptions about F^. The F2 
values are the most precise available from the ZEUS collaboration in the kinematic region 
studied. The extraction of Fl and F2 was based on the reduced double-differential cross 

^ Based on the NNLO calculations by Moch, Vermaseren and Vogt [46,47]. 

3 The conventions used for the CTEQ6.6, ZEUS-,JETS and MSTW08 NLO curves are not the same, for 
example, Fl in CTEQ6.6 is calculated to 0{as) whereas Fl in the ZEUS- JETS and MSTW08 fits are 
calculated to 0{a^) . This accounts for most of the differences. 
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sections, a{x,Q'^), which were measured for 0.09 < y < 0.78 and 20 < < 130 GeV^ 
using data collected at ^/s = 318, 251 and 225 GeV. In addition, and the ratio, R, 
have been extracted as a function of Q^. An overall value oi R = O.I8I005 was extracted 
for the entire kinematic region studied. A wide range of theoretical predictions agree with 
the measured F^. The measurements provide strong evidence of a non-zero F^. 
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-0.6 


0.3 


-0.9 


-1.0 


-0.5 


60 


8.57 


10-" 


0.69 


1.504 


1.9 


3.1 


2.2 


-0.9 


1.0 


-1.3 


-0.2 


0.1 


0.4 


-0.5 


-0.4 


-0.5 


60 


7.89 


10-" 


0.75 


1.586 


2.1 


3.6 


2.2 


-0.9 


-1.7 


-1.6 


-0.6 


-0.4 


-0.3 


-0.5 


-0.9 


-0.6 



Continued on Next Pago. . . 
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X 


y 


-e+p 


&stat 


^sys 










^Ef, 


SelD 


Sdx 


Sdy 


<5mvd 


<5cTD 


<5diff 


(GeV2) 








HER 


(%) 


(%) 




(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


80 


6.07 


10-3 


0.13 


0.884 


0.8 


1.6 




0.6 


0.0 


0.9 


-0.7 


0.0 


0.1 


0.3 


-0.7 


0.5 


-0.1 


80 


3.59 


10-3 


0.22 


1.071 


0.8 


1.7 




0.7 


0.0 


1.5 


-0.2 


0.0 


-0.3 


-0.1 


-0.1 


0.1 


-0.2 


80 


2.54 


10-3 


0.31 


1.204 


1.0 


1.2 




0.8 


0.0 


-0.2 


-0.1 


-0.3 


-0.1 


-0.1 


-0.6 


-0.2 


-0.3 


80 


1.97 


10-3 


0.40 


1.273 


1.2 


1.2 




1.0 


-0.1 


-0.4 


-0.1 


0.0 


0.0 


0.1 


0.2 


0.1 


-0.4 


80 


1.64 


10-3 


0.48 


1.358 


1.4 


1.6 




1.2 


-0.1 


-0.5 


-0.3 


0.1 


-0.1 


-0.3 


0.4 


0.5 


-0.4 


80 


1.41 


10-3 


0.56 


1.463 


1.5 


1.9 




1.4 


-0.2 


0.5 


-0.3 


0.2 


-0.6 


-0.2 


0.3 


0.5 


-0.5 


80 


1.25 


10-3 


0.63 


1.517 


1.9 


2.1 




1.8 


-0.1 


-0.6 


-0.5 


-0.3 


0.4 


-0.1 


-0.3 


-0.2 


-0.5 


80 


1.14 


10-3 


0.69 


1.419 


2.2 


3.2 




2.5 


-0.8 


0.1 


-0.8 


0.3 


-0.1 


-0.5 


0.4 


1.3 


-0.4 


80 


1.05 


10-3 


0.75 


1.436 


2.5 


4.6 




3.2 


-1.2 


1.9 


-1.9 


0.6 


-0.4 


-0.2 


0.9 


-0.9 


-0.5 


110 


8.34 


10-3 


0.13 


0.837 


0.9 


1.3 




0.8 


0.0 


0.7 


-0.7 


0.0 


-0.2 


0.2 


-0.1 


0.2 


-0.1 


110 


4.93 


10-3 


0.22 


1.005 


1.0 


1.7 




0.9 


0.0 


1.4 


-0.2 


-0.1 


-0.2 


-0.2 


0.1 


0.1 


-0.2 


110 


3.50 


10-3 


0.31 


1.157 


1.2 


1.5 




1.0 


0.0 


-0.5 


-0.1 


0.1 


0.3 


0.2 


-0.7 


-0.6 


-0.3 


110 


2.71 


10-3 


0.40 


1.240 


1.4 


1.5 




1.2 


-0.1 


-0.7 


0.1 


0.0 


-0.3 


0.1 


0.2 


0.2 


-0.4 


110 


2.26 


10-3 


0.48 


1.300 


1.6 


1.9 




1.5 


-0.2 


-0.4 


0.0 


-0.1 


-0.4 


0.3 


-0.8 


-0.3 


-0.4 


110 


1.94 


10-3 


0.56 


1.420 


1.8 


2.0 




1.7 


-0.2 


0.6 


-0.2 


0.1 


-0.7 


-0.2 


0.0 


-0.5 


-0.4 


110 


1.72 


10-3 


0.63 


1.409 


2.3 


3.2 




2.4 


-0.4 


-1.6 


-0.3 


0.5 


0.5 


0.1 


0.3 


-0.7 


-0.4 


110 


1.57 


10-3 


0.69 


1.584 


2.6 


3.5 




2.8 


-0.7 


-0.8 


-0.6 


0.7 


-0.2 


-0.3 


0.4 


1.4 


-0.4 


110 


1.45 


10-3 


0.75 


1.717 


4.1 


5.5 




4.6 


-0.9 


1.5 


0.9 


1.2 


0.8 


-0.1 


-0.4 


-1.7 


-0.4 



Table 1: The reduced cross section, a, for the reaction e~^p —>■ e^X at ^/s = 318 GeV . 
The first three columns contain the bin centres in Q"^, x and y, the next three con- 
tain the measured cross section, the statistical uncertainty and the total systematic 
uncertainty, respectively. The final ten columns list the uncorrelated, 6unc o-nd the 
bin-to-bin correlated uncertainties from each systematic source, 5^^, Se^, Se^, ^eiD, 
Sdx, ^dy, ^MVD; ^CTD, ^diflF- For details, see Section\^ A further ±2.7% systematic 
normalisation uncertainty is not included, of which ±2.5% is correlated between the 
running periods and ±1.1% is uncorrelated. 
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Q2 




X 


y 


+ 


Sstat 


^sys 


^unc 


S-yp 


5b. 






^dx 


Sdy 


<5mvd 


<5cTD 


<5diff 


(GeV^) 








MER 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


{%) 


(%) 


(%) 


24 


2.91 


10-3 


0.13 


0.955 


2.4 


3.3 


1.3 


0.0 


0.9 


-0.7 


0.0 


1.1 


0.5 


-1.6 


2.0 


-0.2 


24 


1.72 


10-3 


0.22 


1.105 


2.1 


3.0 


1.1 


0.0 


1.6 


-0.3 


0.2 


1.4 


-0.3 


-1.4 


1.1 


-0.3 


24 


1.22 


10-3 


0.31 


1.153 


2.1 


2.1 


1.0 


0.0 


-0.6 


-0.2 


0.0 


-0.3 


-0.6 


-0.6 


1.4 


-0.3 


24 


9.47 


10-" 


0.40 


1.267 


2.1 


2.3 


1.0 


-0.1 


0.1 


-0.2 


0.5 


1.1 


0.2 


0.8 


1.3 


-0.4 


24 


7.89 


lo-'' 


0.48 


1.328 


2.2 


2.6 


1.2 


-0.3 


-0.7 


-0.4 


-0.4 


-0.3 


-0.2 


0.7 


1.9 


-0.4 


24 


6.76 


lo-'' 


0.56 


1.319 


2.4 


5.0 


1.3 


-0.5 


0.7 


-0.8 


0.2 


-0.5 


0.4 


2.8 


3.7 


-0.4 


24 


6.01 


lo-'' 


0.63 


1.334 


3.0 


3.8 


1.6 


-0.7 


-0.3 


-1.4 


-0.9 


1.2 


0.6 


2.0 


1.6 


-0.5 


24 


5.49 


10-"' 


0.69 


1.432 


3.2 


3.9 


1.8 


-1.0 


-0.5 


-2.0 


-0.2 


-0.3 


0.6 


0.8 


2.2 


-0.6 


24 


5.05 


10-" 


0.75 


1.389 


3.9 


5.0 


2.5 


-2.0 


-0.8 


-3.5 


0.9 


0.2 


0.4 


-0.4 


0.8 


-0.6 


32 


3.88 


10-3 


0.13 


0.909 


1.7 


2.1 


0.8 


0.0 


1.1 


-0.8 


0.0 


-0.2 


-0.1 


-0.3 


1.3 


-0.2 


32 


2.29 


10-3 


0.22 


1.074 


1.7 


2.4 


0.8 


0.0 


1.6 


-0.3 


0.1 


-0.8 


0.2 


-0.6 


1.1 


-0.2 


32 


1.63 


10-3 


0.31 


1.153 


1.9 


2.9 


0.9 


-0.1 


-0.3 


-0.1 


0.0 


0.8 


0.9 


1.0 


2.3 


-0.4 


32 


1.26 


10-3 


0.40 


1.226 


2.1 


2.2 


1.0 


-0.1 


-0.3 


-0.2 


0.4 


0.1 


0.2 


0.6 


1.6 


-0.4 


32 


1.05 


10-3 


0.48 


1.326 


2.3 


1.9 


1.2 


-0.2 


-0.6 


-0.3 


-0.1 


0.2 


0.2 


-0.2 


1.2 


-0.4 


32 


9.01 


10-4 


0.56 


1.270 


2.7 


2.5 


1.4 


-0.5 


0.6 


-0.8 


-0.8 


-0.8 


-0.9 


-0.5 


0.7 


-0.4 


32 


8.01 


lo-'' 


0.63 


1.381 


3.4 


2.9 


1.9 


-0.5 


0.6 


-1.0 


0.6 


-1.6 


-0.2 


0.4 


0.5 


-0.6 


32 


7.32 


lo-'' 


0.69 


1.314 


3.9 


3.7 


2.2 


-0.9 


0.5 


-1.7 


0.4 


0.1 


-0.8 


-1.7 


1.0 


-0.5 


32 


6.73 


10-"' 


0.75 


1.355 


4.5 


5.0 


2.7 


-1.6 


1.6 


-2.8 


-1.1 


0.9 


-0.3 


-1.6 


0.3 


-0.6 


45 


5.46 


10-3 


0.13 


0.890 


1.5 


2.4 


0.7 


0.0 


1.2 


-0.7 


0.0 


-0.3 


0.3 


0.3 


1.7 


-0.1 


45 


3.23 


10-3 


0.22 


1.037 


1.7 


2.5 


0.8 


0.0 


1.2 


-0.2 


0.0 


-0.4 


0.2 


0.9 


1.8 


-0.2 


45 


2.29 


10-3 


0.31 


1.126 


1.9 


1.8 


0.9 


0.0 


0.2 


-0.1 


-0.2 


-1.0 


-0.5 


-0.5 


0.8 


-0.3 


45 


1.77 


10-3 


0.40 


1.171 


2.3 


2.5 


1.1 


-0.1 


-0.4 


0.1 


0.3 


-1.5 


-0.1 


0.6 


1.4 


-0.4 


45 


1.48 


10-3 


0.48 


1.270 


2.7 


2.3 


1.3 


-0.2 


-0.5 


-0.3 


0.3 


0.7 


0.5 


0.4 


1.4 


-0.5 


45 


1.27 


10-3 


0.56 


1.323 


3.0 


2.1 


1.5 


-0.3 


-0.5 


-0.6 


0.1 


0.4 


-0.8 


-0.3 


0.6 


-0.5 


45 


1.13 


10-3 


0.63 


1.385 


3.8 


3.3 


2.2 


-0.6 


-0.8 


-1.1 


1.8 


0.4 


0.5 


0.2 


0.3 


-0.4 


45 


1.03 


10-3 


0.69 


1.443 


4.2 


3.7 


2.4 


-0.8 


-0.6 


-1.3 


0.2 


1.1 


0.4 


1.2 


1.6 


-0.6 


45 


9.47 


10-4 


0.75 


1.400 


4.7 


3.8 


2.6 


-1.0 


0.7 


-1.9 


-1.0 


0.6 


0.8 


0.6 


0.5 


-0.4 


60 


7.28 


10-3 


0.13 


0.809 


1.8 


1.9 


0.8 


0.0 


0.9 


-0.7 


0.0 


-0.2 


-0.1 


0.4 


1.2 


-0.1 


60 


4.30 


10-3 


0.22 


0.970 


2.0 


2.2 


0.9 


0.0 


1.4 


-0.2 


0.0 


-0.2 


0.5 


0.7 


1.1 


-0.3 


60 


3.05 


10-3 


0.31 


1.123 


2.3 


1.8 


1.1 


0.0 


-0.6 


-0.1 


0.0 


-0.1 


-0.4 


0.8 


0.9 


-0.3 


60 


2.37 


10-3 


0.40 


1.183 


2.7 


1.6 


1.3 


0.0 


-0.2 


-0.1 


0.0 


0.7 


-0.3 


-0.2 


0.3 


-0.4 


60 


1.97 


10-3 


0.48 


1.135 


3.3 


2.2 


1.5 


-0.1 


0.6 


-0.5 


0.3 


0.7 


-0.9 


-0.2 


-0.5 


-0.4 


60 


1.69 


10-3 


0.56 


1.277 


3.6 


2.9 


1.8 


-0.3 


-1.0 


-0.3 


-0.1 


1.0 


0.7 


1.4 


0.7 


-0.4 


60 


1.50 


10-3 


0.63 


1.417 


4.4 


3.0 


2.3 


-0.3 


1.3 


-0.4 


0.2 


0.6 


0.4 


0.8 


-0.3 


-0.5 


60 


1.37 


10-3 


0.69 


1.300 


5.2 


4.4 


3.1 


-1.0 


1.8 


-0.7 


0.8 


1.5 


-1.2 


0.5 


-0.4 


-0.4 


60 


1.26 


10-3 


0.75 


1.446 


5.6 


4.8 


3.1 


-0.9 


-1.4 


-2.1 


-0.4 


-1.8 


-0.7 


0.4 


-1.1 


-0.5 



Continued on Next Page. . . 
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X 


y 


-e+p 


&stat 


Ssys 


Sunc 








SelD 




Sdy 


<5mvd 


<5cTD 


5difl 


(GeV2) 








MER 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


80 


9.71 


■ 10-3 


0.13 


0.751 


2.1 


1.7 


1.0 


0.0 


0.8 


-0.7 


0.0 


0.7 


0.2 


-0.3 


0.6 


-0.1 


80 


5.74 


• 10-3 


0.22 


0.947 


2.3 


1.9 


1.1 


0.0 


1.3 


-0.1 


0.1 


-0.7 


-0.2 


-0.5 


0.3 


-0.2 


80 


4.07 


• 10-3 


0.31 


1.067 


2.7 


1.4 


1.3 


0.0 


0.4 


0.0 


0.0 


0.4 


-0.1 


-0.1 


0.3 


-0.3 


80 


3.16 


• 10-3 


0.40 


1.106 


3.2 


2.2 


1.5 


0.0 


-1.4 


-0.1 


0.3 


-0.3 


0.5 


0.1 


0.0 


-0.3 


80 


2.63 


• 10-3 


0.48 


1.170 


3.7 


2.5 


1.8 


-0.2 


-0.7 


-0.3 


-0.5 


-1.0 


0.3 


-0.6 


-0.8 


-0.4 


80 


2.25 


• 10-3 


0.56 


1.217 


4.2 


3.4 


2.2 


-0.4 


-0.7 


-0.6 


0.7 


-1.1 


-0.4 


-1.4 


-1.3 


-0.3 


80 


2.00 


• 10-3 


0.63 


1.246 


5.5 


3.9 


2.9 


-0.5 


0.7 


-1.3 


0.6 


-0.2 


-0.4 


1.3 


1.2 


-0.4 


80 


1.83 


■ 10-3 


0.69 


1.274 


6.1 


4.5 


3.7 


-1.0 


0.3 


-0.4 


0.8 


-1.2 


0.7 


-0.5 


1.5 


-0.5 


80 


1.68 


■ 10-3 


0.75 


1.461 


6.3 


5.3 


3.6 


-0.8 


-2.4 


-1.8 


-1.2 


0.4 


-0.1 


-0.6 


-1.6 


-0.6 


110 


1.33 


• 10-2 


0.13 


0.730 


2.4 


2.1 


1.1 


0.0 


0.8 


-0.6 


0.1 


-0.7 


0.2 


1.0 


0.7 


-0.1 


110 


7.89 


■ 10-3 


0.22 


0.835 


2.8 


2.3 


1.3 


0.0 


1.4 


-0.2 


0.1 


1.1 


0.6 


-0.2 


-0.4 


-0.1 


110 


5.60 


■ 10-3 


0.31 


0.971 


3.2 


1.7 


1.5 


0.0 


0.4 


-0.1 


-0.1 


0.2 


0.4 


0.1 


-0.3 


-0.2 


110 


4.34 


• 10-3 


0.40 


1.078 


3.8 


2.5 


1.8 


0.0 


-0.6 


-0.1 


0.6 


-0.6 


-0.7 


-0.5 


-1.1 


-0.3 


110 


3.62 


• 10-3 


0.48 


1.116 


4.4 


3.5 


2.3 


-0.2 


-2.3 


0.0 


0.0 


-0.3 


-0.3 


-1.1 


0.8 


-0.4 


110 


3.10 


• 10-3 


0.56 


1.192 


5.0 


4.5 


2.7 


-0.5 


-1.0 


-0.4 


-1.5 


1.3 


-0.4 


0.6 


2.7 


-0.5 


110 


2.75 


■ 10-3 


0.63 


1.127 


6.6 


4.0 


3.4 


-0.4 


1.0 


-0.5 


0.4 


-0.9 


0.7 


1.2 


0.3 


-0.5 


110 


2.51 


• 10-3 


0.69 


1.174 


7.7 


6.4 


4.1 


-0.7 


-2.7 


-1.1 


0.4 


-0.6 


-0.4 


-2.4 


3.0 


-0.2 


110 


2.31 


• 10-3 


0.75 


1.354 


10.9 


9.1 


6.2 


-1.0 


4.4 


-1.6 


0.4 


3.3 


0.3 


-1.6 


-2.7 


-0.6 



Table 2: The reduced cross section, a, for the reaction e~^p — ^ e~^X at ^/s = 251 GeV . 
Further details as described in caption of TableUi 
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Q2 




X 


y 


+ 


Sstat 


^sys 


^unc 




<5i5e 








Sdy 


<5mvd 


i5cTD 




(GeV^) 








LER 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


{%) 


(%) 


(%) 


24 


3.64 


■ 10-3 


0.13 


0.864 


1.8 


2.7 


0.9 


0.0 


1.2 


-0.8 


0.0 


0.7 


0.7 


-1.2 


1.4 


-0.1 


24 


2.15 


• 10-3 


0.22 


1.043 


1.6 


2.8 


0.8 


0.0 


1.7 


-0.2 


0.0 


-0.6 


0.3 


-0.6 


1.8 


-0.2 


24 


1.53 


• 10-3 


0.31 


1.136 


1.5 


2.0 


0.7 


0.0 


-0.4 


-0.2 


0.0 


-0.4 


-0.3 


-1.3 


1.2 


-0.4 


24 


1.18 


■ 10-3 


0.40 


1.184 


1.6 


2.0 


0.7 


-0.1 


-0.7 


-0.3 


0.2 


-1.0 


0.3 


-0.6 


1.2 


-0.4 


24 


9.86 


• 10-4 


0.48 


1.195 


1.7 


2.7 


0.8 


-0.2 


-0.2 


-0.4 


0.1 


-0.3 


-0.2 


1.4 


2.1 


-0.4 


24 


8.45 


• lo-'' 


0.56 


1.260 


1.7 


2.0 


0.9 


-0.3 


0.7 


-0.7 


0.4 


-0.8 


-0.2 


-0.8 


0.7 


-0.5 


24 


7.51 


■ lo-'' 


0.63 


1.256 


2.2 


3.1 


1.3 


-0.8 


0.3 


-1.5 
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Table 3: The reduced cross section, a, for the reaction e~^p e^X at ^/s = 225 GeV . 
Further details as described in caption of TableUl 



19 







X 








^2 


^2 


(GeV2) 
















24 


6.67 


■ 10 


-4 




n si+014 


1 4n'?+0 025 

l.'±UO_g 


1 404+0029 
l.'iU'i^g 022 


24 


8.16 


• 10 


-4 


0.25t°;i^ 


27+0-1'* 


1 Q40+0029 


1 o/i 1+0.030 
-'^■'3^^-0.028 


24 


1.08 


• 10 


-3 


0.461^22 


0.4710-?° 


1 260+°-°2^ 

±.Z,UU_g g2g 


1 ORI +0.035 
"— 0.020 


32 


8.89 


• 10 


~4 


0.18t°}l 


0-20l°:J^7 


1 QQ7+0.023 
l.OJ ( _g.023 


1 OQO+0.028 

1.030_g oig 


32 


1.09 


• 10 


-3 


0.26t°;i^ 


27+°-" 

U-^/ -0.10 


1 09^+0.021 
^■'3^^-0.021 


1 09/I+O.O23 

-"^■■^^^-o.oio 


32 


1.43 


• 10 


-3 


^•^ ' -0.19 


0.28l°:?2 


1 ooQ+0.023 

l.ZZ»_g Q23 


1 OOQ+0027 

l.ZZJ_g 016 


45 


1.25 


• 10 


-3 


14+0" 

'-'•-'^^-0.11 


15+0-12 


1 QO/1 +0.024 

-L-OZ4:_g Q24 


1 00/1+0.029 
^■'3^^-0.018 


45 


1.53 


• 10 


~3 


-0 


o-ool°:J° 


-1 900+0.022 

l.ZOO_g Q22 


1 246+°-°2'^ 
^■^^"-0.010 


45 


2.02 


• 10 


-3 


' -0.19 


0.381°:?^ 


1 170+0.022 

l.i < o_g.022 


1 170+0.028 

l.i 1 0_g 016 


60 


1.67 


• 10 


-3 


^•-'^'^-0.12 


17+013 

^•-■^ ' -0.09 


1 O9R+0.029 

-L-OZO„g 029 


1 OOK+0 034 

-L-OZO„g 021 


60 


2.04 


• 10 


-3 


n 19+015 


o.2iirii 


1 9-1 -1 +0.025 
l.Zll_g 025 


1 01 1 +0.030 
l.Zil_g 017 


60 


2.69 


• 10 


-3 


0.271°;- 


0.28l°:?l 


1 14t;+0 024 
i.l^O_g 024 


1 145+°-°^° 

1.140_g 015 


80 


2.22 


• 10 


-3 


0.181°;!^ 


19+0-15 

'-'•-'^■^-0.09 


1 25fi+°-°23 

l.ZdU_g 029 


1 256+°-°3'* 

l.ZOD_g 022 


80 


2.72 


• 10 


-3 


0.381°;}^ 


0.391°:}° 


1 91 O+0.027 

l.Z10_g 027 


1 O1/I+0.032 

l.Zl'i_g 022 


80 


3.59 


• 10 


-3 


"•-^^-0.22 


13+019 


1 nS9+0 022 
-L-UOZ g 022 


1 n«l+0 029 

i.U8i_g 010 


110 


3.06 


• 10 


-3 


33+015 


S4+0-18 

'-'•'5^-0.12 


-, 9t-4+0.033 
l.Z01_g 033 


1 of;c;+0 039 

1-Zdd_g 02g 


110 


3.74 


• 10 


-3 


0.311^;}^ 




1 1 Qfi+0.029 
J^-loU_g.029 


1 1 OQ+0.035 
l-lOO_g 023 


110 


4.93 


• 10 


-3 


O.Oll°;i 


^•"-3-0.03 


1 noo+0 026 

l.UZZ_g 026 


1 noo+0.037 
l.UZZ_g 006 



Table 4: T/ie measured values of Fl and F2 at 18 separate {x,Q'^) points. The 
quoted uncertainties include both the statistical and systematic sources, whereas a 
±2.5% correlated normalisation uncertainty is not included on the and F2 val- 
ues. The (1) superscript indicates an unconstrained fit whereas the (2) superscript 
indicates that the constraints F2 > and < F^ < F2 were enforced by the prior. 
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Table 5: The single values of Fl extracted in each bin. The quantities are 
quoted such that Q'^/x was constant for each value, which corresponds to y = 0.71 
for ^/s = 225 GeV. The quoted uncertainties include both the statistical and sys- 
tematic sources, although a ±2.5% normalisation uncertainty is not included. The 
(1) superscript indicates an unconstrained fit whereas the (2) superscript indicates 
that the constraints F2 > and < F^ < F2 were enforced in the prior. 
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Table 6: The single values of R extracted in each bin. Other details as in the 
caption to TablelE, although no additional ±2.5% normalisation uncertainty need 
be included. 
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Figure 1: Detector-level distributions of the energy, E'^, and polar angle, 9^, of the 
scattered electron candidates within the HER, MER and LER data sets compared to 
the combined MC predictions (MC DIS+BG). The background only MC is labelled 
MC BG. The vertical dashed-line represents the E'^ cut. The 9^ distributions are 
shown for E'^ > 6 GeV. 
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Figure 2: The reduced cross sections at 6 values of Q'^ as a function of x for the 
three centre- of-mass energies. The points represent the ZEUS data from the HER 
(n), MER (a) and LER (•), respectively. The solid lines represent the predicted 
reduced cross sections, using the ZEUS-JETS PDFs. The dashed lines represent 
the predicted reduced cross sections when Fl is set to zero. The points and lines 
are shifted by c (see top right) for clarity. The inner error bars represent the sta- 
tistical uncertainty. The outer error bars represent the statistical plus systematic 
uncertainties added in quadrature. A further ±2.7% systematic normalisation un- 
certainty is not included. 



23 




Figure 3: and F2 at 6 values of as a function of x. The points represent the 
ZEUS data for Fi (•) and F2 (k), respectively. The error bars on the data represent 
the combined statistical and systematic uncertainties. The error bars on F2 are 
smaller than the symbols. A further ±2.5% correlated normalisation uncertainty is 
not included. The DGLAP-predictions for F^ and F2 using the ZEUS-JETS PDFs 
are also shown. The bands indicate the uncertainty in the predictions. 
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Figure 4: Values of (a) Fl and (b) R as a function of . The error bars 
on the data represent the combined statistical and systematic uncertainties. A 
further ±2.5% correlated normalisation uncertainty is not included. The shaded 
band labelled ZEUS average represents the 68% probability interval for the overall R. 
The lines represent various model predictions (see text for details). 
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